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Abstract 

We have measured the branching fraction of the inclusive radiative B meson decay 
B^Xgj to be 

BriB^Xsj) = (3.36 ± 0.53(stat) ± 0.42(sys)l|jJ^(th)) x 10"^. 

The result is based on a sample of 6.07 x 10^ BB events collected at the T{4S) 
resonance with the Belle detector at the KEKB asymmetric e+e~ storage ring. 



PACS: 13.40.Hq, 14.40.Nd, 14.65.Fy 
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In the framework of the Standard Model (SM), the inclusive decay B^Xg^y, 
where Xg is a hadronic recoil system containing an s quark, proceeds primarily 
through the electroweak penguin diagrams of the radiative b^sj transition. 
These diagrams are of interest because they have loops where non-SM par- 
ticles, such as charged Higgs bosons or SUSY particles, can contribute and 
potentially produce deviations from SM expectations. An SM branching frac- 
tion calculation for B—^Xg^ that includes next-to-lcading order QCD correc- 
tions has an expected precision of 10% [1]. Experimental measurements of 
the branching fraction at this level of precision are useful for identifying or 
limiting non-SM theories [2]. 

A semi- inclusive analysis is used to reconstruct the B^Xg'y decay from a 
primary photon, a kaon and multiple pions. The monochromatic photon energy 
from the two-body decay b^sj is smeared by gluon emission, the Fermi motion 
of the b quark inside the B meson, and the B meson momentum in the T(45') 
center-of-mass (CM) frame. 

The B^Xgj decay is studied with the Belle detector at the KEKB asym- 
metric e+e" storage ring [3]. The dataset consists of a sample of 6.07 x 10^ 
BB events corresponding to an integrated luminosity of 5.8 fb~^ taken at the 
T(4S') resonance {^/s = 10.58 GeV), and an off-resonance background sample 
of 0.6 fb"^ taken 60 MeV below T(45). The b eam energies at the T(4S') reso- 
nance are 3.5 GeV for positrons and 8.0 GcV for electrons. Quantities used in 
this analysis are defined in the CM frame, in which the beam energy -Ebeam is 
defined as -\/i/2. Inclusion of charge conjugated modes are implied throughout 
the text. 

A full description of the Belle detector is given in [4]; here we briefly describe 
the components relevant for this analysis. Charged tracks are reconstructed 
inside a 1.5 T magnetic field induced by a super-conducting solenoid magnet, 
by means of a 50 layer central drift chamber (CDC) that covers the labora- 
tory polar angle 17° < 6iab < 150°. Tracks are flttcd through the CDC and a 
three layer double sided silicon vertex detector (SVD), with a transverse mo- 
mentum resolution of {(Tpjptf = (0.0019pt)2 + (0.0034)^ {p^ in GeV/c). For 
particle identification, we combine information from three detectors. The high 
momentum range, typically from 1 to 3.5 GeV/c, is covered by a silica aero- 
gel Cherenkov counter (ACC) system, providing threshold-type kaon-to-pion 
separation. The momentum range below 1.5 GeV/c is covered by a time-of- 
flight (TOF) counter system, and the very low and high momentum ranges 
are accessible with dE/dx information from the CDC. Photons are detected 
by an electromagnetic calorimeter (ECL) located between the particle identifi- 
cation devices and the solenoid coil. The entire tracking acceptance is covered 
with 8736 Csl(Tl) crystals that are typically 5.5 x 5.5 cm^ in cross-section 
at the front surface and 16.2 Xq in depth. The photon energy resolution is 
{aE/Ef = 0.0132 + (0.0007/E)2 + (0.008/^1/4)2 -^^ q^^), which is better 
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than 2% for the primary photons used in this analysis. The Belle detector 
is modelled with the GEANT program [5] to simulate the detector response. 
Both the data and the Monte Carlo (MC) events [6] are reconstructed with 
the same program. 

Hadronic events are selected with an efficiency of 99% for BB final states. 
Events with at least three tracks from the interaction region and visible energy 
greater than0.2y^ are selected. The QED and r+r^ events that remain are 
removed using ECL information and event topology [7]. 

Photon candidates are selected from ECL clusters of 5 x 5 crystals. Each pho- 
ton candidate is required to have a laboratory energy greater than 20 MeV 
with no associated charged track, and shower shape consistent with an elec- 
tromagnetic shower. We accept events in which the most energetic photon is 
in the barrel region (33° < 9iab < 132°). We require that 95% of its energy 
is concentrated in the central 3x3 crystals. The efficiency of these criteria 
for the primary photon is (80 ±4)%. In order to reject 7r° and t] mesons, we 
combine this primary photon candidate with any other photon (for r] we only 
use photons with > 200 MeV) and reject the event if the invariant mass 
of the two photons is within ±18 MeV/c^ of M^o or ±32 MeV/c^ of M^. 

We reconstruct the recoil system Xg in 16 difi'erent final states of one charged 
kaon or Kg plus one to four pions which may include one 7r°. We combine the 
primary photon with every Xg combination and calculate the opening angle 
dxs-y between and px^ , the energy difference AE — E^ + Ex^ — Eheam ^-nd 
the beam constrained mass M^c — \/ E^eam ~ IPt + Px^ | ^ ■ Then, we require 
Ox^-y > 167°, Mbc > 5.2 GeV/c^ and -0.15 < AE/GeV < 0.1. 

Charged kaon candidates are selected using a kaon-to-pion likelihood ratio. 
For every charged track, hkelihoods from ACC, TOP and dE/dx are calcu- 
lated for pion and kaon hypotheses. A combined likehhood is constructed for 
each hypothesis, and a tight cut is applied on the likelihood ratio. Kaons are 
identified with a typical efficiency of 65% with 2% pion fake rate. All charged 
tracks that are not identified as kaons are considered to be pions. 

The neutral kaon candidates are reconstructed in the — > tt+tt" mode 
from two oppositely charged tracks. The Kg momentum is recalculated with a 
vertex constrained fit. The candidates are then required to have an invariant 
mass within ±8 MeV/c^ (~ 2(t) of M^o and a vertex that is displaced from 
the interaction point in a direction consistent with the Kg momentum. The 
Kg reconstruction efficiency is (65 ± 4)%. 

The 7r° 77 candidates are reconstructed from two photons that have an 
invariant mass within ±16 MeV/c^ (~ 3(j) of M^o. Then, the tt*^ momentum 
is recalculated with a mass constrained fit. The 7r° reconstruction efficiency is 
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(60 ± 4)%. 



When multiple candidates are found in an event, the best candidate is se- 
lected as follows: If there is at least one charged kaon or pion that forms the 
Xs vertex with the constraint from the run-by-run determined profile of the 
B meson decay vertices [8], the candidate with the largest vertex confidence 
level is chosen. The candidate with the largest Oxs'y is chosen in the follow- 
ing cases: (1) there is no charged kaon or pion to form the vertex; (2) there 
remains an ambiguity whether to add a 7r°; or (3) there are several 7r° and 
Kg candidates to choose from. After selecting the best candidate, we require 
Mx, < 2.05 GeV/c^. 

The largest background source is from continuum qq production where the 
photon originates from initial state radiation or from high momentum neutral 
hadron decays such as 7r° — > 77 with one undetected photon. In order to 
reduce and estimate the contribution of the qq background, we introduce a 
new variable that exploits the topological difference between the jet-like qq 
events and the spherical BB events. We define the following variables in the 
B meson rest frame, 

^ ^ J2i \Pi\\p^\Pi{cosei^) ^ ^ Eij \Pi\\pj\Pi{coseij) 

' Ei\Pi\\p^\ ' ' Ei,j\Pi\\Pj\ 

where Pi is a Lcgcndrc polynomial of order /, and i, j run over all the photons 
and charged tracks that are not used to form the B candidate. Then, we 
combine six of them into a Fisher discriminant 

F^Yl ^i^i + E An, 
1=2,4 1=1,2,3,4 

where the six coefficients ai and f3i are optimized to maximize the discrimi- 
nation between signal and background as shown in Fig. 1. The terms Ri and 
R3 are excluded from the Fisher discriminant, because either of these terms 
is found to have a correlation with M^c- We call F the Super Fox-Wolfram 
(SFW) variable, since the terms are combined in such a way as to enhance 
the discriminating power of the original Fox- Wolfram moments [9] . We select 
events with F > 0.1, which is a compromise between the statistical significance 
of the signal and the size of the systematic error due to background subtrac- 
tion. In order to estimate the qq background contribution, we use the SFW 
sideband region of F < —1.5, where the signal fraction is only (0.7 ± 0.2)%. 
This translates into a contribution of 0.7 ± 0.2 events in the signal region. 

In addition to the qq background, there is a contribution from 6— >c decays, 
in which B^D'^*^p modes are dominant. This background, which does not 
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contribute to the SFW sideband, is estimated by MC to be 9.1 ± 1.8 events 
and is subtracted from the data sample. 

The signal yield is obtained from the M^c distribution (Fig. 2). In Fig. 3 we 
show that the ratio between the number of events in the signal region and 
in the sideband of the SFW variable is independent of the beam constrained 
mass Mbc according to a MC calculation assuming contimmm qq production. 
This observation is consistent with the background estimated using the off- 
resonance data. Thus, the SFW sideband data can be used to model the 
background shape of the M^c spectrum in the absence of a substantial sample 
of off-resonance data. We fit the Af^c distribution with a background shape 
taken from the SFW sideband data and a signal shape obtained from MC 
simulation. 

We observe 222 events in the signal region M^c > 5.27 GeV/c^, of which 
115.5 ± 7.7 is the estimated background contribution, and is consistent with 
the MC expectation. An excess of 106.5 ± 16.8 events is clearly seen in the 
distribution shown in Fig. 2 at the B meson mass with a mass resolution 
of 4.4 MeV/c^ and a signal-to-background ratio of 0.9. 

To obtain the background subtracted recoil mass spectrum shown in Fig. 4, we 
use the shape of the Mx^ spectrum determined from MC, whose contribution 
is normalized to the estimated background of the M^c distribution. We adopt 
this method because the ratio between the number of signal events and SFW 
sideband events has a finite slope as a function of Mx, , especially in the region 
> 2.05 GeV, and thus the SFW sideband does not properly represent the 
background sample for the Mx^ spectrum [10]. Nevertheless, the slope in the 
Mxs < 2.05 GcV region is small and compatible with a fiat distribution within 
the statistical fluctuation, and the background subtraction method using the 
SFW sideband sample leads to consistent results for the Mx^ spectrum and 
the S— >Xs7 branching fraction. 

The signal MC sample is generated as a mixture of the exclusive A'* (892)7 
mode and an inclusive B^Xg'f contribution for Mx^ > 1.15 GeV/c^ in order 
to separate the K*{892) signal from higher resonances starting with Xi(1270). 
The recoil system Xg is modelled as an equal mixture of sd and su states [6]. 
To simulate the inclusive Mx^ spectrum, we adopt the model by Kagan and 
Neubert [11] with the b quark pole mass parameter rrib = 4.75 GeV/c^. 

The fraction of K*(892)'y decays in the B^Xg^ transition, rK*-y, is deter- 
mined from the Mx, spectrum with a consideration of the migration effects 
on the reconstructed Mx^ due to incorrect particle assignments to the Xg 
final state. The Mx^ spectrum is divided into two regions: the A'* (892) re- 
gion below 1.15 GeV/c^ and the continuum region between 1.15 GeV/c^ and 
2.05 GeV/c^. The MC sample is analyzed to determine the number of events 
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migrating between the regions. The MC results are used to unfold the migra- 
tion effects and to generate a set with an improved value, rx*'y = (f f .3±3.5)%. 
Using the new MC sample, the procedures of M^c fitting, Mx^ background 
subtraction, and rx'-y unfolding are iterated. It is checked that rx*7 is stable 
under the iteration. The size of the migration effects and the value of rx*7 
depend on the assumed value: for = 4.75 GeV/c^, 77% of the events 
reconstructed below 2.05 GeV/c^ are genuine and 23% are from the region 
above 2.05 GeV/c^. We note that due to a non-negligible amount of migration 
between the regions below and above 2.05 GeV/c^, this boundary does not 
correspond to a sharp value in Mx^ or E^. 

According to the authors of ref. [11], the K*{%%2) contribution can be approx- 
imated by integrating the Mx^ spectrum up to a threshold mass above the 
-ft'*(892) mass, though their model does not explicitly take resonance contri- 
butions into account. This integral is sensitive to the parameter rrib. To match 
the measured ii'*(892) fraction we varied mf, and the threshold mass and ob- 
served that rrib > 4.9 GeV/c^ is disfavored. We use nib = (4.75 ±0.10) GeV/c^ 
as the parameter region to test the model dependence of the branching frac- 
tion measurement [12]. We assume a somewhat larger error on m^, than that 
used for the semileptonic B meson decays in the combined LEP, SLD and 
GDF evaluation [13] based on the recent theoretical work of [14], in which 
mb{l GeV) = (4.58 ±0.06) GeV/c^ is quoted. We note that the corresponding 
pole mass is nib — 4.70 GeV/c^. Other theoretical input parameters have little 
impact on the shape of the Mx^ spectrum. 

The event selection efficiency is determined from a MC sample that is cali- 
brated with high statistics control data samples for all final state particles. 
The reconstruction efficiency for the primary photon is obtained using pho- 
tons with energies between 2 and 3 GeV from the processes e+e~^e+e~7 and 
fj^jj. The charged tracking efficiency for momenta relevant to this analysis 
is obtained by comparing the high momentum r] yields for the rj^-K^-K'Ti^ 
and ?7^77 decay modes. The kaon selection efficiency is determined from 
a sample of (p^K^K^ and Dg^cpir decays. The pion rejection efficiency is 
also determined from the rj^n'^n'n^ sample. The Kg and 7r° efficiencies 
are tested with a sample of K*{S92) decays into K~^7i~, K'^ir'^, Kg7r~^ and 
Kgii^. Combining all 16 different channels, the selection efficiency becomes 
(12.8±1.3)%. The Mx,, SFW and Mbc cuts reduce the efficiency by (57±2)%, 
(50 ± 3)% and (70 ± 1)%, respectively The SFW cut and the n^/rj veto effi- 
ciencies are measured by applying the cuts on a sample of exclusively recon- 
structed >D°7r~, D'^^K~7r'^ decays, in which the first pion is assumed 
to be massless in order to mimic the primary photon for the cuts. The size 
of the systematic error from the efficiency determination is dominated by the 
statistical error of the control sample. The final event reconstruction efficiency 
is (2.58 ± 0.29)%. 
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The total systematic error is quoted as a quadratic sum of the systematic 
errors on the efficiency, the best candidate selection procedure, signal MC 
statistics, SFW sideband uniformity and the number of BB events. The sensi- 
tivity of the cuts to the experimental error on the beam constrained mass M{,c 
is tested by changing the MC resolution function, and the effect is found to be 
negligible. The efficiency errors on the 16 different final states are calculated 
individually and are then combined. The best candidate selection is tested 
with a B~^D^7:~ sample, which is analyzed using the inclusive reconstruc- 
tion method of this analysis. The effect of the SFW sideband uniformity is 
tested by adding a slope to the background distribution within the tolerance 
level of Fig. 3. The contributions to the systematic error on the branching 
fraction are given in Table 1. 

The fitted signal yield is extrapolated over the entire Mx^ region using the 
model of ref. [11] with = 4.75 GeV/c^ and corrected for the final states 
that are not reconstructed. The analysis is repeated for rrih = 4.65 GeV/c^ 
and 4.85 GeV/c^ to test the model uncertainty. The signal efficiency is found 
to be 15% lower and 16% higher, respectively. This error and the error of r^*^, 
which typically corresponds to 5.0% in the branching fraction calculation, are 
combined into the theoretical model error. 

Finally, we obtain the branching fraction of B^Xg^, 

Br(B^Xs-f) = (3.36 ± 0.53(stat) ± 0.42(sys)t^J^(th)) x 10"^ 

where the first error is statistical, the second is systematic and the third is the 
theoretical model error. Using the values of the measured B^Xg'j branching 
fraction and the fraction of the i^* (892)7, we find the B i^* (892)7 branch- 
ing fraction to be (3.8 ±0.9) x 10^^ (statistical error only), which is consistent 
with the B — > i^* (892)7 branching fraction measurements from CLEO [15] 
and a separate Belle analysis [16]. In Fig. 5 we show the measured photon 
spectrum in which the background is subtracted and the efficiency loss due 
to the cut on Mx^ is corrected by using MC input. This spectrum may be 
compared with model calculations. 

To summarize, we have measured the inclusive branching fraction of B^Xg^ 
decay with the Belle detector. The result is consistent with the SM prediction 
of Br{B^Xs-i) = (3.28 ± 0.33) x 10"^ [1]. The result is also consistent with 
previous measurements by the CLEO [17] and ALEPH [18] experiments. 
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Table 1 

The contributions to the systematic error on the branching fraction. 



Photon reconstruction 


±5.3% 


Charged track reconstruction 


±4.7% 


Charged kaon selection 


±1.8% 


Charged pion selection 


±1.1% 


Kg reconstruction 


±1.2% 


7r° reconstruction 


±3.1% 


SFW efficiency and Tr^/ry veto 


±6.8% 


Best candidate selection 


±2.4% 


SFW sideband uniformity 


±4.7% 


Signal MC statistics 


±3.0% 


Number of BB events 


±2.3% 



Total systematic error ±12.4% 




Fig. 1. The distribution of the SFW variable described in the text. The off-resonance 
background data (solid circles) and the qq MC expectation (open histogram) are 
compared with the signal distribution of B^Dir data (open circles) and signal MC 
(hatched histogram). 
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2 

Beam constrained mass (GeV/c ) 

Fig. 2. The beam constrained mass (M(,c) distribution (a) compared with the to- 
tal background qq and 6— >c (open histogram); (b) after background subtraction 
compared with the signal MC expectation (hatched histogram). 




Beam constrained mass (GeV/c^) 



Fig. 3. The ratio of the background events in the SFW signal region to the sideband 
region as a function of M^c from the qq MC (solid circles) and off-resonance data 
(open circles). 
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Recoil mass (GeV/c ) 



Fig. 4. The recoil mass {MxJ distribution after background subtraction compared 
with the signal MC expectation (hatched histogram). 
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Fig. 5. The photon energy spectrum, background subtracted and corrected for the 
cut-off on Mxg- The data points are compared with signal MC expectations for 
three different values of m^. 
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